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During the evolution of the universe there are at least two epochs during which electromagnetic 
waves cannot scan the universe's internal structure neither bring information to outside observers. 
The first epoch is when photons are in local thermodynamic equilibrium with other particles, and 
the second is when photon scattering by charged particles is strong. One can call these two periods 
of cosmological time as standard unlighted epochs. After the last scattering surface, photons become 
relic photons and turn into a source of information about the universe. Unlighted cosmic epochs 
can also appear when one considers non-minimal theories, i.e., theories in which the electromagnetic 
field is coupled in an intricate way with the cosmological gravitational field. By considering a 
cosmological model where the dark sector, i.e., the dark energy and dark matter, self-interacts 
via an Archimedean-type force, and taking into account a non-minimal coupling theory for the 
electromagnetic field, we discuss the appearance of unlighted epochs. In the framework of our non- 
minimal theory, a three-parameter non-minimal Einstein-Maxwell model, the curvature coupling can 
be formulated in terms of an effective refraction index n(t). Then, taking advantage of a well-known 
classical analogy, namely, in a medium with n 2 < electromagnetic waves do not propagate and their 
group velocity, i.e., energy transfer velocity, has zero value at the boundary of the corresponding 
zone, one can search for the unlighted epochs arising in the interacting dark fluid cosmological 
model. We study here, both analytically and numerically, cosmological models admitting unlighted 
epochs. 
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I. INTRODUCTION 
A. The cosmology, two-component dark fluid and Archimedean-type interaction 

In order to examine important features happening within the universe, such as light propagation, one has to set up 
from the start a model for the dynamics of the universe as a whole, smoothing out all irregularities. General relativity 
tells that both geometry and matter are important in the dynamics of the universe and, moreover, through Einstein's 
equations one finds that geometry guides the matter and matter changes the geometry. 

In relation to geometry, we now have a good idea of the spacetime geometry of the universe. It is governed by 
a cosmological time t, which is the proper time of the fundamental particles, usually considered as galaxies, of the 
substratum. Slices of this time yield, by the use of the cosmological principle, a homogeneous and isotropic spatial 
geometry, so that the metric is the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric. From observations the 
universe is expanding, which is then taken into account by a single function, the cosmological scale factor a(t), which 
once known, also yields the rate of expansion, i.e., the Hubble function H{t) = |, and the rate of acceleration, i.e., 
the acceleration parameter —q(t) = -^fjj- Moreover, the spatial slices are flat, or almost flat, simplifying the problem 
even further (see, e.g., [1]). 

In relation to matter, we also have now a good idea of the matter content of the universe. There are three main 
components, namely, dark energy which accounts for 70% of the matter, dark matter which accounts for around 
25%, and the rest, which includes baryonic matter, radiation, and other forms such as black holes, which accounts 
for 5%. Dark energy and dark mater are essential building blocks in any cosmological model. Their precise nature, 
both of the dark energy and dark matter, is unknown, but this does not preclude a good understanding of the gross 
features of the dynamics of the universe. Dark energy 0-U is essential in the explanation of the observed accelerated 
expansion of the universe @ . Dark matter provides an excellent explanation for the flat velocity curves seen in 
the outskirts of a galaxy and for the dispersion velocities of galaxies in clusters of galaxies @, [1] • Given there are 
these two important matter components that essentially drive the dynamics, it is worth contemplate them as a single 
whole. There are various manners in which this could be performed. One intriguing possibility is considering both 
components as manifestations of a single dark fluid j9l4l3j|. Other possibility is to postulate an interaction between 
them so that the two components although really distinct have to be treated in a broader unified scheme, see e.g., 

MM- 

In [16l Il7j this latter possibility of a unified interaction scheme between dark energy and dark matter has been 
followed. The background for the interaction itself was given by relativistic hydrodynamics theory for the dark energy 
component and a relativistic kinetics framework for the dark matter component. Baryonic matter is negligible in 
the context of cosmological dynamics and so it has been left out in the scheme. The interaction itself between 
both components has been provided by an Archimedean-type four-force, a direct generalization of the Archimedean 
buoyancy law, in which the three-force is proportional to the gradient of the pressure. The dark energy pressure 
is of the same order as the dark energy energy density and for this reason, the influence of the Archimedean-type 
force on the dark matter component can be important. Now, the dark energy pressure can be negative, in which 
case the Archimedean-type force can be negative, instead of positive as in the usual buoyancy force case. In the 
case the pressure changes sign as the universe evolves, there appear several different stages in its acceleration, in 
which the acceleration itself can change sign. It was further shown in [ll| [lj| that the Archimedean-type force is 
able to distribute between the two fluid components the corresponding energy content of the universe, which in turn 
guides the whole evolution of the universe, in its one or several stages of acceleration. Multistage universes have also 
appeared in [ill ]. 

B. Non-minimal coupling 

The way in which the electromagnetic field couples to gravity is an open question. Usually it is assumed a minimal 
coupling where simple flat spacetime derivatives involving the field are replaced by covariant derivatives but this 
might not be so. For instance, one can replace the flat spacetime derivatives by covariant derivatives along with terms 
involving the curvature tensor and its contractions, yielding a non-minimal coupled theory. In a gravitational strong 
regime, like in the vicinity of black holes or in a cosmological context, these additional terms, if present, can be felt 
and are important (see, e.g., original papers and reviews [l9|-|4(|). 

A simple and most fruitful non-minimal theory, a non-minimal Einstein-Maxwell theory, possesses an action with 
a Lagrangian which includes a linear combination of three cross-invariants, namely, qiRF mn F mn , qiR %k Fi m Fk n g mn 
and q3R lkmn FikF mn , where q\, qi and q^ are free parameters. The scalars appearing in the Lagrangian are linear in 
the Riemann tensor R l kmn , Ricci tensor Run, and Ricci scalar R, and are quadratic in the Maxwell tensor Fa-. An 
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interesting fact about this non-minimal theory is that the corresponding non-minimal electrodynamic equations have 
the same form as the equations for the electrodynamics of anisotropic inhomogeneous continuous media (see, e.g., 
[4ll - l43l |). This circumstance makes it possible to consider non- minimal analogs of well-known phenomena in classical 
electrodynamics of continuous media, such as birefringence induced by curvature (30j , anomalous behavior of electric 
and magnetic fields non- minimally coupled to the gravitational radiation curvature induced Cherenkov effects 
, and optical activity in vacuum [46] . Of course, the variety of non-minimal effects depends on the spacetime of 
the model and its symmetries. 

Since electromagnetic waves, from radio waves, to light, to gamma rays, are a very important source of information 
about the universe, it is worth to study the changes that might operate in a universe with a non-minimal coupling 
between the electromagnetic and gravitational fields. 

C. Non-minimal electromagnetic wave propagation in an expanding universe with a cosmic dark fluid with 

Archimedean-type interaction and unlighted epochs 

1. Prologue 

In a generic gravitational background field, electromagnetic waves non-minimally coupled propagate with a velocity 
which differs from the velocity of light in vacuum. These waves, influenced by tidal interactions induced by curvature, 
do not travel along null geodesies of the background spacetime. Thus, non-minimal coupling may produce significant 
changes in the propagation of electromagnetic waves. It is therefore worth to study this coupling in a cosmological 
context. 

In cosmology, in particular, when one deals with a spatially homogeneous isotropic expansion for the universe, one 
is faced with electromagnetic effects, possibly non-minimal, of two types: first, the phase and group wave velocities 
depend on the cosmological time, and second, the amplitude and energy density of the electromagnetic waves decrease 
with time. Different aspects of these two phenomena have already been studied within a non-minimal setting. For 
instance, in (38j the deviation of the photon velocity from the velocity of light in vacuum was estimated in the 
framework of Drummond and Hathrell's approach |2l| . 

2. Motivation 

We are interested here in developing the study of the propagation of electromagnetic waves non-minimally coupled 
to the gravitational field in an isotropically expand ing universe. Since cosmological models with dark fluid interaction 
are of great interest, such as those studied in [l6l . Il7j . there are now new motives to study anew the propagation of 
these electromagnetic waves non-minimally coupled to the gravitational field in an isotropically expanding universe 
governed by a dark fluid. We mention three such motives. 

The first motive is related with the synergy between non-minimal electromagnetic wave propagation and a cosmology 
in which a dark fluid (i.e., a fluid with two components, namely, dark energy plus dark matter) plays a major role 
[l6l . In our non-minimal theory there are three coupling constants q\ , qi and c/3 . These non-minimal constants 
can be determined directly from a more fundamental theory, as was done in the case of quantum electrodynamics 
vacuum polarization effects in a curved background by Drummond and Hathrell [2l|, or they can be considered as 
phenomenological inputs. Of course, the results for q\, q2 and 93, in plj should be thought as providing a portion 
of the whole values to the three constants, since other effects can provide curvature induced interactions between the 
fields, and all the effects should be summed together. Our idea here is that a dark fluid can also support non-minimal 
interactions. Thus, we consider the parameters q\, qi and (73 to be phenomenologically introduced and do trust that 
a mechanism of non-minimal coupling between electromagnetism and gravity via a dark fluid can be clarified in the 
future. In other words, if a dark fluid, as a medium with unusual properties, can act as a mediator of non-minimal 
coupling between photons and gravitons, we expect that the coupling parameters q\ 1 qi and q^ will be estimated 
theoretically and established from observations. 

The second motive is related to the existence of unlighted epochs in our models. The non-minimal coupling of the 
electromagnetic waves with the gravitational field can be codified in terms of an effective time-dependent refraction 
index n(t), where t denotes cosmological time. From this, we will show that for specific values of the non-minimal 
parameters q±, q2 and 53, n 2 (t) can be negative during some finite time interval, and can either vanish or take infinite 
values at some transition points t( s ), i.e., points where the sign of n 2 (t) changes, the subscript s, labeling the several 
possible transition points s = 1 , 2, ... and denoting sign change. This means that the universe can pass through epochs 
when electromagnetic waves cannot propagate as their phase velocity is a pure imaginary quantity during this period 
of time. Taking into account that in order to read the history of the universe one should reconstruct the sequence of 
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the events by using the whole spectrum of electromagnetic radiation, we coin such epochs as unlighted epochs, since 
portraits of the universe of those epochs cannot be available. Moreover, the equations for trapped surfaces in the 
universe have here the form t = £( s ), and they are time-like in contrast to the standard space-like ones that appear 
in black hole formation. The use of the term unlighted epoch is then well fit, since it distinguishes clearly the two 
situations, namely unlighted epochs in the universe versus trapped surfaces in black holes. 

The third motive is determined by the necessity to divide properly the history of the universe into epochs according 
to various physical scenaria. Since the transition points in the universe history are fixed by the critical temperatures 
Tr c \ of some basic physical processes, one should link the temperature evolution T(t) and the cosmological time t by 
using the cosmological scale factor a(t). When the cosmological time scale is defined upon using relic cosmic microwave 
background photons traveling along null geodesic lines, the corresponding law is very simple, T(t)a(t) — T(t )a(t ), 
where t is some reference time, like now, say. On the other hand, when we take into account a non-minimal coupling 
between photons and gravitons, the corresponding link between the temperature T{t) and the scale factor a(t) is much 
more sophisticated. Indeed, it is determined not only by a(t), but also by its first and second derivatives, a(t) and 
a(t), respectively. In other words, the temperature T(t) is tied to the scale factor a(t), the Hubble function H(t) = ^ 
and the acceleration parameter — q(t) = -^p- This fact, of course, should induce novel and interesting features. 

3. This work 

In [l6l . [l7j a cosmological model based on a dark fluid made of two components, the dark energy, considered as the 
fluid substratum, and the dark matter, described in the framework of kinetic theory, interacting via an Archimedean- 
type force was considered. The evolution models for the universe were classified in terms of the admissible transition 
points, i.e., the points in which an accelerated expansion of the universe is changed to a decelerated expansion and 
vice versa. This classification is based on the analysis of the function acceleration parameter — q(t) = -^jp-- Now, we 
find here that the function n (t), describing the behavior of the effective refraction index, also includes the functions 
H(t) and —q(t). Thus, in order to analyze non-minimal light propagation, it is useful to use n 2 (t) instead of —q(t). In 
this work we study, analytically and numerically, the propagation of electromagnetic waves non-minimally coupled, 
along with their phase and group velocities, in six cosmological models based on a dark fluid which self interacts 
through an Archimedean- type force [l|| H3] ■ This work is thus an interesting sequel of pil [l?} • 

4- The organization of the paper 

The paper is organized as follows. In Section|TT]we introduce the cosmological model. In subsection lll Al we give the 
cosmological model itself and recall the basic formulas of the Archimedean-type model, which are necessary for further 
numerical analysis. In subsection III Bl we give the basic set up for non- minimal interaction and light propagation. In 
Section III CI using the master equations of non-minimal electrodynamics we reconstruct the effective dielectric and 
magnetic permeabilities, the effective refraction index n(t), the effective (optical) metric in terms of the scale factor 
a(t), the Hubble function H(t), acceleration parameter —q(t) and the three non-minimal coupling phenomcnological 
constants q 1: q 2 and q^. We also discuss the refraction index n(t), focus on the derivation of the expression for the 
group velocity of the waves (i.e., the velocity of energy transfer) using the analogy with electrodynamics of continua, 
and define unlighted epochs. In section UTTI we make an analytical study of the unlighted epochs. In subsection IIII Al 
using the Kohlrausch stretched exponential functions we give the cosmological models fit for the study, in subsection 
IIII B I we provide the choice for the non-minimal models, and in IIII CI we give some analytical cosmological examples 
describing unlighted epochs. In Section HVl we consider the results of numerical analysis for the refraction index n(t), 
phase and group velocities and effective lengths of the photon trips for six basic Archimedean-type submodels. In 
subsection IIV Al we discuss the cosmology, in subsection IIV Bl we provide the choice for the non-minimal models, and 
in subsection IIV CI the models are analyzed in detail, including the submodels of perpetually accelerated universe, 
periodic and quasiperiodic universes, and various submodels with one, two and three transition points. We also study 
the photon pathlength, the true duration of epochs and the universe lifetime for these models. In Section [V] we draw 
conclusions. In the Appendix we analyze a number of one-parameter examples for which the coupling constants are 
linked by two relations motivated physically and geometrically. 
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II. THE COSMOLOGICAL MODEL 
A. Cosmological context 

1. The Cosmology 



We start from the action functional 



S = I d^xj—g ( ^ + L matter ) , ( 1 ) 



2 k 

where R is the Ricci scalar, g is the determinant of the four-dimensional spacetime metric gik, L mattcr is the matter 
Lagrangian which we assume includes a cosmological constant A, and k=^^- with G being Newton's constant and 
c the velocity of the light. Using the standard variation procedure with respect to the metric g lk one can obtain 
Einstein's equations. Latin indices run from to 3. In the above we assume that £ matter , anc j so the corresponding 
energy-momentum tensor, T J ™ attcr , comes from both the dark energy and dark matter, which together make up for 
most of the total energy of the universe. Thus, the master equations for the gravity field are assumed to be the usual 
Einstein's equations, 

Ri k -\Rg ik = KT^, (2) 

where Rik is the Ricci tensor, given by the contraction of the Riemann curvature tensor R\ mk . 

We use the spatially homogeneous flat metric of Friedmann-Lemaitre-Robertson- Walker (FLRW) type given by 

ds 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) , (3) 

where we chose units with c=l. The energy momentum tensor is described by the functions p(t) and LT(t) , and E(t) 
and P(t). The functions p(t) and IT(i) describe the energy density and pressure of the dark energy, respectively. 
The functions E(t) and P(t) are the energy density and pressure of the dark matter, respectively. The cosmological 
constant A is incorporated into the dark energy state functions, i.e., the term is included into p and II (la ]. With 
these assumptions the equations coming out of Eq. ([2]) for the gravitational field read 



d\ 2 8ttG 



a / 3 



(P + E), (4) 



v a J \a J 3 
with the dot denoting a derivative with respect to time. 



[(p + E) + 3(U + P)}, (5) 



2. Two-component dark fluid: Archimedean-type interaction between dark energy and dark matter 

In order to describe the evolution of several quantities, in particular of the refraction index n(t), as this is the 
quantity we are interested, we need of the Hubble function H(t), defined as 

H(t) EE ° , (6) 

a 



and the acceleration parameter —q(t), 



Sometimes it is useful to swap —q(t) for H, given by, 

H(t) = -H 2 (t)[l + q(t)]. (8) 
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Now, we obtain these functions using the Archimedean-type interaction between the dark energy and dark matter 
model 0,113 • The function H(t) can be found from Eq. © and Einstein's equation (QJ, yielding 

H\t) = ^(p + E). (9) 
The acceleration parameter —q(t) can be found from Eq. ([7} and Einstein's equation ([5]), yielding 

(10) 

As for the function H(t), one uses Eq. © and Eqs. ([§])- (|T0|) to hnd 

H(t) = -4irG (p + E + II + P) . (11) 

Thus, we need to find the following state functions: the energy density of the dark energy p(t) and its pressure 

the energy density of the dark matter E(t) and its pressure P(t). The search scheme for these quantities is the 

following. First of all introducing a new convenient variable 

"ITS' (12) 
a(t ) 

and using the auxiliary formulas 

d „. . d f o<*o) dx , „. 

xH{x) t -t = -—- (13) 

we find, through some manipulation of Einstein's equations a key equation for Tl(x) 

^x 2 U"(x)+xW(x) (4£+a) +3{l+a)U+3p ^J(x) , (14) 

a second order differential equation, linear in the first and second derivatives, II' = -4-11 and II", and nonlinear in the 
function Tl(x). The parameters £, a and po are the coupling constants involved into the assumed linear inhomogeneous 
equation of state, namely, 

p(x) =p + aTL(x)+Zx-^IL(x). (15) 

This equation links the pressure to the energy density of the dark energy (see [l6j for details). The source-term 
J(x) = J(x, II— 11(1), IT) in the right-hand side of (|14[) is given by the integral 



where we used the definitions 



( ir (a) — ^ — v^mm » (16) 



F (a)(*) = 4exp{2V (a) [n(l)-n( : r)]}, (17) 



_ jV(a)m (a) A (a) _ m (a) 

^(a) = TFT\ ^ ' A ( a ) = 1 ^ ' 



/>oo 

if l/ (A( a )) = / cosh i/z exp [— A( a ) cosh z] 
Jo 



(19) 



Here V( a ) are the constants of the Archimedean- type coupling (see [16| for details), K v {\^) are the modified Bessel 
functions of second kind of order v (we are interested in the v = 2 case), 11(1) = H(t ) is the initial value of the dark 
energy pressure, /cb is the Boltzmann constant, and T( a ), rnu), and Nu>, are the partial temperature, the mass, and 
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the number of particles per unit volume of the sort (a), respectively. When the quantity n(ic) is found, other state 
functions can be calculated using the integrals 

E ( x ) = E "nr / y 2 dyJi+y 2 F (a) (x) e - A <*> (20) 



P(x) = y f" e -^V^ (21) 

Thus we have formulated our integration scheme. From here, we want to work out novel features that appear in such 
a cosmological context. 



B. Non-minimal coupling: set up 

Now we want to study light propagation envisaged as a perturbation in the background cosmological manifold. 
This means that the electromagnetic field propagates in a given geometry, without modifying the geometry itself. 
The electromagnetic field is a test field. This is well justified since it is known that the energy density of radiation is 
thoroughly negligible in relation to the dark energy, dark matter and baryonic matter. We also assume that the field 
is non-minimally coupled and write the electromagnetic action functional as, 

^electromag = 1 j {F mn F mn + TZ ikmn F ik F mn ) . (22) 

Using the standard variation procedure with respect to the electromagnetic potential four-vector A4 one can obtain 
the non-minimally extended Maxwell equations (see, e.g., 43], see also 34]). Here = diA k ~d k Ai is the Maxwell 
tensor. Then, the non-minimal Maxwell equations have the following form 

V fe iT fe = 0, V k F* ik = 0, (23) 

where the excitation tensor H lk is linked with the Maxwell tensor F mn by a linear constitutive equation 

H lk = C ikmn F mn , (24) 

with the linear response tensor C lkmn given by 

The dual tensor F* lk = |e ' mn F mn is defined in a standard way through the Levi-Civita tensor e zkmn . The tensor 
■j^ikmn j s non-minimal susceptibility tensor decomposed as 

n i kmn _ ^( 9 ™y"-/" s ta )+|(ff™(, fc "-r/ ra +i?' !n 9 m -]? tm j' n )+ g 3ff' : ~ , (26) 

where qi , qi , and qz are the phenomcnological parameters of the non- minimal coupling, and R lkmn ; R lm t and R, are 
the Riemann tensor, the Ricci tensor, and the Ricci scalar of the background cosmological manifold, respectively. It 
is appropriate here to recall that the contracted susceptibility tensor 7V m satisfies 

K im = g kn K lkmn = l -Rg im {Z qi + <? 2 ) + R m (q2 + 5a) , (27) 

and so vanishes in a generic curved spacetime when the non-minimal coupling parameters are linked through two 
relations 3qi+q2 = and 92 + 93 = 0. Analogously, the scalar 1Z given by 

K = g mi9kn l¥ kmn = g m K m - R{Qq x + 3q 2 + q 3 ) , (28) 

has zero value in a generic curved spacetime when 6gi+3q2+<?3=0. 
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C. Non-minimal light propagation in a cosmological context: Optical metric and refraction index, phase 

and group velocities, and unlighted epochs 

1. Optical metric and refraction index 

The optical metric - Based on the metric ([3]) and on the symmetries of the spacetime one obtains that the non- vanishing 
components of the non-minimal susceptibility tensor lZ lk mn , see Eq. ([26)) . have the following form 

2 

a , 



-K U lt = -K 2t lt = -n U lt = {iq 1 +2q 2 + q3 )- + {^i+q2) (-) , (29) 



ft 12 = -ft is = -ft 23=(3«l+«2)-+(3«l+2«! + <? 3 ) 



(30) 



where the indices 1,2,3 correspond to x,y,z, respectively. In addition, due to the spacetime isotropy the linear 
response tensor C lkmn can be rewritten in the standard multiplicative form 



C 



■/ ki. 



i 



(gi m 9 



kn in km\ 

y* y* ) 



where 



n\t)8\8? 



a 2 (t) 



(d 1 d 1 +d 2 d 2 +d 3 d : 



3 ; i 



is the so-called associated metric (see [13]), and the function n 2 (t) is defined as 

n 2 (t) = e(t)n(t), 

with e{t) and [i(t) given by 

s(t) ee 1 + 2K u lt , = i + 2TZ 12 12 . 



(31) 



(32) 



(33) 



(34) 



Taking into account that in a cosmological context t he g lobal velocity four- vector is U l =61, one can see that the 
associated metric (|3"2"j) is in fact an optical metric (48l. |49| given by 



= g im +[n 2 (t)-l]U i U n 



Its inverse is then 



with 



9kr, 



9kr. 



n 2 {t) 



- 1 



UkU„ 



Am * _ ri 
9* 9km — °k 



(35) 



(36) 



(37) 



holding. The quantity n(t) is then interpreted as an effective refraction index due to its correspondence, in conjunction 
with the optical metric gl m , to the effective refraction in the geometrical optics approximation framework. 

Let us see this in detail. In the geometrical optics approximation, the electromagnetic potential Ak and the field 
strength Fki can be put as follows 



A k = A k e™ 



F H =i 



p k Ai -piA k 



(38) 



for some amplitude A k and phase function ty. In the leading order approximation the Maxwell equations (|23[) reduce 
to C* kmn Pk Pm A n = 0, i.e., 



[" im kn n in „km 

I y* y* 



9T9*\ PkPmA n =0. 



(39) 
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Using the optical metric property gl m U m — n 2 U\ the Landau gauge U m A m — 0, and the requirement that the 
frequency u) = U m p m is non- vanishing, one can show that Eq. (|39[) leads to the dispersion equation 

9 k r Pk Pm = . (40) 

Thus, the propagation of photons non-minimally coupled to the gr avity field is equivalent to their motion along a null 
geodesic line in an effective spacetime with metric (|35|) (see, e.g., [48l l49l| for details). In our setting, the refraction 
index n(t) determines the effective phase velocity of light in our FLRW-type spacetime, namely, V p h = ^fey- Thus, 

let us consider in more detail the properties of the quantity n 2 (t). 

The expression for the refraction index - The square of the effective refraction index can be obtained from Eqs. (|33[1 - 
(|3~4"|) and is given by 

n 2 (t) = l-2(3 gl +2 g2 +g 3 )|-2(3 gl + g2 )(f) 2 

1 - 2(3^+52)1 - 2(3^+2(72+53) (f) 2 ' 

Notably, it can be put as a function of the Hubble function H(t) and the acceleration parameter —q(t), defined in 
©-([7]), as well as of two effective non-minimal coupling constants, Q\ and Q 2 given by, 

Qi = -2(3gi+2g 2 +g 3 ) , Q2 = -2(3 qi +q 2 ) . (42) 
As a function of these quantities Eq. (|4T|) can be rewritten as 

2 _ l + [Q2-Qiq(t)}H 2 (t) 

U[) l + {Ql-Q2q(t)}H 2 (t)- l4dj 

For the cosmological models with a de Sitter-type final stage (i.e., — q(t —> oo) — > 1), Eq. (|43|) yields n 2 —> 1. The signs 
and the values of the parameters Q\ and Q 2 are not yet established, but one can discuss several phenomenological 
possibilities using geometric analogies and physical motivations, see Appendix. 



2. Phase velocity and group velocity (or energy transfer velocity) 
The phase velocity of an electromagnetic wave is defined as (c = 1) 

v °^l = W)- (44> 

This definition is standard and follows directly from the dispersion relation [50j . 

On the other hand, the definition of group velocity, or energy transfer velocity, V* , is connected with the definition 
of the electromagnetic field energy flux, and this problem requires a preliminary discussion. First of all, one should 
state that formally the propagation of electromagnetic waves takes place in isotropic backgrounds with a refraction 
index equal to one, but under the influence of a non-minimal coupling the interaction can be reformulated in terms of 
an effective refraction index n(t), which depends on time through the Riemann tensor components. Since, as we have 
seen, the non-minimal coupling of photons to gravity is equivalent to the consideration of some isotropic dispersive 
medium, it seems reasonable to take the relevant stress-energy tensor of the electromagnetic field for the definition 
of the energy transfer velocity V* . Now, there is a number of definitions for the electromagnetic stress-energy tensor 
in a medium, the best known are the ones by Minkowski [5l| , Abraham [52], HH , Grot-Eringen-Israel-Maugin (53 - 
[57j . Hehl-Obukhov 58], and de Groot-Suttorp (5t|. The energy flux four- vectors in these definitions differ from one 
another. Thus for us a choice between the several definitions is considered as an ansatz. 

We follow the Hehl-Obukhov definition [58| . according to which the stress-energy tensor formally is as in vacuum, 
^Tfc leCtI ° maS ' i- e -' T° A l° ctIomas = \gikF mn F mn — F im F k "\ and thus is symmetric, traceless and does not depend on the 
macroscopic velocity of the medium. Nevertheless, the ponderomotive force 

pi ^ ^rpclcctvom^gkl pi ^^,(^/^ m pkm^ (45) 

is non- vanishing. In contrast to the vacuum case, this stress-energy tensor is not a conserved quantity, its four- 
divergence is equal to the ponderomotive force. According to this definition [581 ] . the energy flux P is 

P = (g ip - [/ » C /P) T clcctromag [/g _ rfPq BpEq , (46) 
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where 



jf pq = e ipqs U. 



Bp 



FLU*, E q = F qj W, 



and the energy-density scalar W is 



W = UPT° q ectmmas U q = -^{E m E m + B m B m ) . 



(47) 



(48) 



When we deal with test electromagnetic waves in an isotropic universe, the necessary quantities can be found as 
follows. The non-minimally extended Maxwell equation can be rewritten now in the form 



V fe 



1 



l 



(F* m U k - F k m W)U r ' 



0, 



(49) 



where fi(t) and n 2 (t) are given by (f3"3")) - (j3"4")) . Let the direction along which the electromagnetic wave propagates be 
x , say. Then Eq. (|49| admits the following solution for the potential four- vector, 



with 



A l {t,x 1 ) = Sf 6(t)cos* , 
* = *(*o) + fell* 1 -/(*)], 



(50) 
(51) 



and where k\ is a constant wave-vector, b(t) the amplitude factor, and f(t) the frequency function. These satisfy the 
following equations, 



4f 
dt 



an 





d 


( an 2 ^ 


db~ 




dt 


A V J 


~dl_ 



(52) 



d_ 

dt 



b 2 an' 



( 1 
dt 



0. 



(53) 



In the geometrical optics approximation, i.e., for short waves, one has (fci — > oo) and (|52]) gives / = i.e., the 
quantity u)(t) = k\f = ^ plays the role of a time dependent frequency, the quantity k(t) = ^ is a wave- vector 
modulus, and the dispersion relation takes the form uj(t) = ^^. In this approximation, Eq. (1531 gives b(t) = b ^J 

From a physical point of view (see, e.g., (50[) it is reasonable to calculate the non- vanishing energy flux four- vector 
component and the energy density scalar averaged over a wave period, in which case one can use (cos 2 "J) = (sin 2 = 
2 and (cos ^ sin^P) = 0, where ( ) denotes average over a period. In the geometrical optics approximation the energy 
flux P and the energy-density scalar W have the form 



1 



I 1 ) = - — (F a F 12 ) = 



and 



(W) 



1 



■{F 2 2 ) 



1 

2a 2 ^ t2/ ' 2a* 

Thus, the physical component of the energy transfer velocity V* is 



V* 



(W) 



b 2 k\ 
2na 5 



b 2 k\{n 2 + 1) 



2n 



+ 1 



(54) 



(55) 



(56) 



The function V*(n(t)) = n 2 " 1 is appropriate for the description of the energy transfer velocity, since depending on 
the refraction index n it vanishes in the limits n — > and n — ¥ oo. In addition, its value is maximum in pure vacuum 
since V* = 1 when n = 1, and generically satisfies the inequality |V*(n(f))| < 1. Below, to simplify terminology, we 
indicate this energy transfer velocity as group velocity V gT , i.e., V gr = V*, and so 



(57) 
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3. Unlighted cosmological epochs: definition 

It is remarkable that the square of the effective refraction index can take negative values, n 2 (t) < 0, for some set of the 
parameters Qi and Q 2 . We call the time intervals for which n 2 (t) < as unlighted epochs, since during these periods 
of time the refraction index is a pure imaginary quantity, and the phase and group velocities of the electromagnetic 
waves are not defined. The function n 2 (t) can change sign at the moments tu\ of the cosmological time when the 
numerator or the denominator in Eq. (|43l) vanish. When the numerator vanishes, one has n(t/ s \) = 0, V^>h (£(«)) = 00 > 
and V SI (tr s \) = 0- When the denominator vanishes, one has n(t( s )) = oo, Vph(i(s)) = 0, and V gI (tr s \) — 0. In both 
cases the unlighted epochs appear and disappear when the group velocity of the electromagnetic waves vanishes, i.e., 
at these points there is no energy transfer. In our terminology the times t( s ) are the unlighted epochs boundary points. 
For this reason, the condition Vgr(ti s )) = sets the criterion for the unlighted epoch appearance or disappearance. 
In other words, the unlighted epochs start and finish, when the associated metric, i.e., the optical metric, given in 
Eq. (j3"2")l becomes singular. 

In our cosmological context we distinguish three types of unlighted epochs: 

1. Unlighted epochs of the first type: These epochs start at tn)=0 with n 2 (0) < and finish with n 2 (i( 2 ))=0- 

2. Unlighted epochs of the second type: These epochs start at tm > with 7i 2 (i( 1 ))=0 and finish at tm > tm with 

n 2 (i (2 ))-0. 

3. Unlighted epochs of the third type: These epochs appear when at least one boundary point has a refraction index 

with an infinite value, i.e., n 2 (t^)—oo or n 2 (i( 2 ))=o°- Of course, in this type both quantities may be infinite. 

Clearly, only one unlighted epoch of the first type can exist, whereas the number of unlighted epochs of the second 
and third types is predetermined by the guiding parameters Qi and Q 2 . 



III. UNLIGHTED COSMOLOGICAL EPOCHS: ANALYTICAL STUDY 



A. The cosmology 

In order to illustrate analytically the physics of unlighted epochs, let us consider a scale factor a(t) given by a 
stretched exponential function, namely, 

a(t) =a exp{(Ht) u }, (58) 

for some constant % and exponent v. This function was introduced by Kohlrausch [6(| in 1854 and now is sys- 
tematically used in various physical and mathematical contexts (see, e.g., papers concerning the applications of the 
Kohlrausch- Williams- Watts (KWW) function [UH^I). When v=l the function ([551) coincides with the standard de 
Sitter exponent. When v=2 one deals with an anti-Gaussian function obtained in [16j as an exact solution of a model 
with Archimedean-type interaction between dark energy and dark matter. This stretched exponent also appeared in 
in the context of generalized Chaplygin gas models. For the function (|58")l one obtains from Eqs. (JB])-©, 

H{t) = V U u t»-\ - q (t) = l + —{Ut)-\ H(t) = v(v-l)H u t v - 2 . (59) 

v 

It is reasonable to assume that v is positive in order to guarantee that a(t — > oo) — > oo and — q(t — > oo) — > 1. Clearly, 
when < v < 1, H is negative, and the Hubble function H(t) vanishes at t — > oo. Our ansatz is that v < 2. Note 
that for v < 2 one has H(t oo) — > 0, and thus, n 2 (t — > oo) — > 1. Let us focus on the properties of the function 
n 2 (t) for a scale factor of the Kohlrausch type Eq. (|58|) . 



2(f] = (Hi) 2 -" + ^n 2 (Q 1 +Q 2 )(my + v{v-l)U 2 Q 1 

{) (my-" + v 2 u 2 (Q 1 +Q 2 ){uty + v { v -i)u 2 Q 2 ■ 1 ' 

When v = 1, there are no unlighted epochs, since n 2 (t) = 1. When v — 2, both the numerator and the denominator 
in (|60p are quadratic functions of i, which means that depending on the values of the parameters Q\ and Q 2 , one can 
find zero, one, two, three or four boundary points i( s ). In other words, the number of unlighted epochs can be zero, 
one or two. These models represent perpetually accelerated universes since— q(t) is non-negative for t > 0. When v 
is a perfect rational v=^^ with mi and m 2 natural numbers and mi < m 2l the numerator and denominator of the 
function (1601) can be rewritten as a polynomial of order 2m 2 — mi. Such a polynomial has 2m 2 — mi roots, part of 
them, say k, can be real and positive defining the corresponding number of unlighted epochs. 
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FIG. 1: A plot of n 2 (t) as a function of time t in the Kohlrausch-type model. Left panel: Contains sketches of unlighted epochs 
of the first and third types. For both types the Kohlrausch parameter is set to v — 1.00001. The non-minimal parameters 
are the following: Q2 = —Qi = 10000 for the model of the first type, and Q2 = —Qi = —10000 for the model of the third 
type. The unlighted epoch of the first type appears in the plot as a dashed curve. It is continuous, it starts with a negative 
value, n 2 (0) = — 1, and has one zero (the point in which the phase velocity is infinite and the group velocity is vanishing). The 
unlighted epoch of the third type appears in the plot as a solid curve. It is discontinuous, the function n 2 (i) has one infinite 
jump (in this point the phase and group velocities are vanishing). Right panel: It is illustrated a unlighted epoch of the second 
type with parameters v = 2/3, Qi — —81/800, Q2 = —729/800. There are two points, in which the refraction index and group 
velocity vanish and the phase velocity becomes infinite. 



B. The choice for the non-minimal parameters 



In order to develop the unlighted epochs in the cosmological models of interest one has to choose the non-minimal 
parameters. Here we give two example which are enough to have a feeling of the physics. In one case we put Q\ = —Q2, 
and define Q = Q\. In the other case we choose Qi = — oo^La and Q2 ~ 2 ' - 



MoW ^ ~~ soon 2 



C. Two cosmological models and their unlighted epochs 

1. Universe with v < 2 and Q\ = — Q2 = Q 
We now put Q\ = —Q2 and define Q = Q\. Then, the square of the refraction index in Eq. (|60p has now the form 

There is now only one time moment t — when n 2 is either equal to zero or infinite. Equivalently, there is only one 

moment when the group velocity takes zero value. When (y— V)Q < 0, the denominator is positive, n 2 (0) = —1 < 0, 

1 

and the numerator vanishes at = ^ \yT-L 2 \(v— 1)Q|] 2 ~" ; clearly, we deal with an unlighted epoch of the first type, 
see the dashed line in Fig. [T] left panel. When (v—l)Q > 0, the numerator is positive, n 2 (0) = — 1 < 0, and the 

denominator vanishes at tm = -h \y{v—l)'H 2 Q\ 2 ~ u ; clearly, we deal with an unlighted epoch of the third type with 
n2 (*(i)) = °°j see the solid line in Fig. [TJleft panel. 
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2. Universe with v = ^ and Qi — — 80 q^ , Q2 = — 8() ' im j 

Using the auxiliary variable 2; = {2THt)i we can rewrite the square of the refraction index in Eq. (I60[) in the 
following form 

2(] = x 2 +AH 2 (Q 1 +Q 2 )x-18H 2 Q 1 
n{ > x 2 +m 2 {Q 1 +Q 2 )x-im 2 Q 2 ' 1 ' 

It is easy to find general conditions, when the denominator has no roots and the numerator has two positive roots. 
For instance, Fig. [T] right panel, illustrates the case with Qi = — 80 q^ and Q 2 = — gjq-^ ■ We deal now with an 

unlighted epoch of the second type, which starts at i(i)=57 [^(3 — V5)] 2 , and finishes at t^ 2 ) — j^ [^j(3 + V^5)] 2 • 

IV. UNLIGHTED COSMOLOGICAL EPOCHS: NUMERICAL STUDY. THE REFRACTION INDEX, 
PHASE AND GROUP VELOCITIES FOR SIX MODELS AND THE EFFECTIVE PHOTON 

PATHLENGTH 

A. The cosmology 

In [l6l . Il7j several cosmological models were presented through numerical results. Here we study non-minimal light 
propagation in six such models. They include, perpetually accelerated universes, periodic universes, universes with 
one transition, two transition and three transition points, and quasiperiodic universes. 

B. The choice for the non-minimal parameters 

To study non-minimal light propagation in the universes indicated above, there is a plethora of possibilities to make 
a choice of the non-minimal parameters. A sample of those possibilities is presented in the Appendix. For our study 
we make one choice, bearing in mind that the other choices will reproduce qualitatively the same results. We choose 
that non-minimal susceptibility scalar 1Z vanishes, 1Z = 0. From Eq. (|28[) this means 6qi+3q 2 +qs = and thus, from 
Eq. P^]) . one has Qi + Q 2 = 0. Putting Q = Q\ = —Q 2 , one obtains in this one-parameter example the following 
expression for n 2 (t), 

a _ l-QH 2 (t) [l + g(t)} _ 1 + QH(t) 
U[t) l + QH 2 (t)[l + q(t)} l-QH(t)- [M) 

Now that we see an explicit expression for n 2 (t) our choice can be motivated as follows. First, it is an example that 
best illustrates the several unlighted epochs. This is because only one function, H, guides the behavior of the effective 
refraction index. Second, when t — > 00, H ~ > and thus n 2 — > 1 providing V p h — > V gr — > 1 in our late-time Universe, 
as it should be. We will consider this one-parameter example as the one which will provide a substratum for our 
numerical analysis. 

C. Six cosmological models, their unlighted epochs, and the photon pathlength 

The results of our numerical calculations for the six basic cosmological models in which an Archimedean-type 
interaction between dark energy and dark matter plays a main role are presented in Figs. [2HZ1 The panels (a) of the 
figures display the plots of H(t). The panels (b) display the plots of the square of the refraction index n 2 (r), where 
r is a logarithmic time defined by r = logx, with x = ^^y- For the calculations we use formula (|63|) putting there 

the values of the guiding parameter Q for which the square of the refraction index n 2 (t) can be negative for some 
period of time. The panels (c) display the plots of the electromagnetic wave phase velocity. Here we use values of Q, 
for which n 2 (t) > and the phase velocity happens to be a pure real function. The panels (d) display the plots of 
the group velocity for the same values of the guiding parameters Q. Additional internal windows clarify fine details 
of the plots. 
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1. Perpetually accelerated universe 



The first class of models is the class of perpetually accelerated universes. This class arises when — q(t) is non- 
negative for t > and so there are no points in which — q(t) could change sign. In this sense the history of such class 
of universes includes only one epoch. 

In this model the dark energy energy density p is always non-negative and the dark energy II is non-positive, see 
[l7j for more details. Although there is only one epoch for these universes, one can divide this epoch into eras. The 
extrema of the functions —q(t), H(t), p(t), and Tl(t) give the boundary points between the eras, see 17] for more 
details. We do not plot —q(t), H(t), p(t), and n(£), see [17[ for such plots. 

We are interested in light propagation and unlighted epochs. So we plot H(t), the square of the refraction index 
n 2 (t), the phase velocity V p h(t), and for the group velocity V gl -(t). Fig. [2] shows that indeed for the particular case 
chosen there are at least six eras which appear visually. 

The first era is an era of superacceleration. During this era the function H(t) grows monotonically and reaches its 
global maximum, as it is shown in the panel (a) of Fig. [2] The same behavior happens for the square of the refraction 
index n 2 , for the phase velocity Vph(i), and for the group velocity V gI (t), as it is shown in panels (b), (c), and (d) of 
Fig. [21 respectively. 

The second, third, and following eras appear one after the other as the parameters —q(t), H(t), p(t), and II(i) 
of the universe relax to a state with asymptotically constant positive values, namely, — goo, H^, p^ and negative 
IIoo. For the parameters that guide the properties of light propagation we see from Fig. [5] that H(t) asymptotically 
vanishes, see panel (a), and the refraction index, the phase velocity and the group velocity tend quasiperiodically 
asymptotically to one, see panel (b), (c), and (d), respectively. 

Other features can be mentioned: (i) There are at least three eras, in which the refraction index exceeds one, n > 1, 




FIG. 2: Perpetually accelerated universe. The plots in panels (a), (b), (c), and (d) are presented for a typical model with the 
following parameters: £=0.35, cr=-0.99 (i.e., 3£+cr ~ 0.06 > 0), V (0) =l, £ (0) =0.0205, A (0) =l, p,=0.333 • 10" 4 , and II'(1)=-1. 
Panel (a): illustrates the behavior of the function Hit), the plot of which has a finite number of visible damped oscillations, 
and tends to zero asymptotically; it relates to the case of the non-negative acceleration parameter —q(t), i.e., there are no 
transition points in this model, there is no partition of the universe history into epochs, but there are a few eras, the start and 
finish of which are marked by the extrema of the function H. Panel (b): The plot of the square of the refraction index n 2 (r) 
is presented in the panel (b) for the parameter Q — 0.7. It contains at least three eras, in which the refraction index exceeds 
one, three eras with n < 1, and n 2 (r) tends asymptotically to one at t — > oo. There is an unlighted epoch of the first kind, 
which is characterized by negative values of the function n 2 (i); it starts at t = and extends to the middle of the first era. At 
the end of this unlighted epoch the refraction index takes zero value, and starting from this point electromagnetic waves can 
propagate and transfer information into the universe. Panel (c): This panel displays the phase velocity of the electromagnetic 
waves as a function of time; the plot relates to a parameter Q = —0.55; this choice guarantees that n 2 (i) is positive, and V p h 
is a real function everywhere. The plot of the phase velocity reflects the history of universe, i.e., it has extrema just at the 
moments when one era is changing into another. Clearly, the waves move more slowly in the early universe, when the curvature 
is large, then there are few eras with oscillations of the phase velocity near the vacuum speed of light, and finally, Vph tends 
to one asymptotically. Panel (d): The last panel shows the behavior of the group velocity of the electromagnetic waves; the 
calculations are made for the same values of the guiding non-minimal parameters as for the phase velocity. Clearly, Vg T does 
not exceed one, it tends to one asymptotically, and an energy transfer takes place slowly in the early universe. 
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and three eras with n < 1. (ii) There is one unlighted epoch of the first type, which is characterized by negative values 
of the function n 2 {t). This unlighted epoch starts at t — with a non-zero value of n 2 , and extends up to the middle 
of the first era. At the end of this unlighted epoch the refraction index takes zero value (the corresponding effective 
phase velocity would be infinite), and starting from this point the electromagnetic waves can propagate and transfer 
information into the universe, (iii) Clearly, the waves move more slowly in the early universe, when the curvature 
is large, (iv) Let us note that the plot of the phase velocity properly reflects the history of the universe, i.e., it has 
extrema just at the moments when one era is changing into another. 

2. Periodic universe 

The second class of models is the class of periodic universes. This class arises when the equation q(t) — has 
an infinite number of roots, and the history of the universe splits into an infinite number of identical epochs with 
accelerated and decelerated expansions. 

The acceleration parameter —q(t), and the Hubble function H(t), oscillate with fixed frequency and amplitude after 
the second transition point [17j . These oscillations are reflected in the oscillations of the energy density and pressure 
of the dark energy. We do not plot —q(t), H(t), p(t), and II(i), see [ItJ for such plots. 

We are interested in light propagation and unlighted epochs. From the plots Fig. [3] we see that H(t) also oscillates, 
which in turn reflects again in the behavior of the refraction index n 2 (t), the phase velocity V p h(t), and the group 
velocity V gY (t). 

Other features can be mentioned: (i) The plot of the function n 2 (t), see panel (b) of Fig. [3l displays an unlighted 
epoch of the first kind, which extends from t=0 into the middle of the first era. 

3. One transition point universe 

The third class of models is the class of universes with one transition point during their evolution. This class arises 
when the acceleration parameter — q(t) is a deformed Heaviside step-function, or a hyperbolic tangent (see [17[). 

The acceleration parameter —q(t) has the property that the change of a deceleration epoch into an acceleration 
epoch takes place only once in a narrow period of time, and the plot looks like a typical plot for a phase transition. 



H(t) n 2 (x) 




FIG. 3: Periodic universe. The plots in panels (a), (b), (c), and (d) are presented for a typical model with the following 
parameters: £=0.1, a=-0. 299999 (i.e., 3£+<r ~ 10" 6 > 0), V (0 )=l, £ (0 )=0.0205, A (0 )=l, p,=0.333- 10" 4 , and n'(l)=0.01. The 
model is characterized by infinite number of transition points, in which the acceleration parameter — q changes sign. Panel (a): 
The function H(t) (for Q — 15) grows monotonically in the early universe and then oscillates near zero value. Panel (b): The 
plot of the function n 2 (t) displays, first, an unlighted epoch of the first type, then grows monotonically up to the end of the 
first era, and then oscillates near the value n = 1. Panel (c): This is a plot of the phase velocity of the electromagnetic waves 
(for Q = —0.25). The curve follows the plot of H, indicating the starting and finishing points of the corresponding epochs in 
the history of the universe. Panel (d): This is a plot of the group velocity. It shows it also grows monotonically during the first 
era, then starts to oscillate so that its maximal value reaches the speed of light in vacuum. 
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The second epoch, characterized by an accelerated expansion, looks like the de Sitter-type stage with p+H = and 
vanishing E and P. The first and second epochs are not divided into eras within this model. We do not plot —q(t), 
H(t), pit), and Il(i), see [13] for such plots. We are interested in light propagation and unlighted epochs. Fig. 0] 
illustrates this class of models. 



4- Two transition points universe 

The fourth class of models is the class of universes with two transition points and one extremal point during their 
evolution. This class arises when the acceleration parameter — q(t) has a double root, i.e., an extremal point, inside 
the second epoch of acceleration. This model is also distinguished by the fact that the dark matter pressure reaches 
a maximum at the end of the first era of the first acceleration epoch. We do not plot —q(t), Hit), pit), and TI(t) , see 
pj} for such plots. 

We are interested in light propagation and unlighted epochs. Fig.[5]illustrates this class of models. The sophisticated 
behavior that appears in this model is reflected explicitly in the form of the plot H(t), see Fig. [SJ In this scenario 
there is an unlighted epoch of the second kind, where the function n 2 (t) vanishes at some t, with t > 0, then takes 
negative values at the end of the first era of the first epoch, and then vanishes again. Thus this unlighted epoch is 
separated by two points, where the refraction index vanishes and the effective phase velocity takes infinite values. 

5. Three transition points universe 

The fifth class of models is the class of universes with two epochs of decelerated expansion and two epochs of 
accelerated evolution. This class arises when the acceleration parameter —q(t) is of the so-called TV-type [17] . This 
model is also distinguished by the fact that the start of the universe's expansion relates to the deceleration epoch, 
which is then replaced by a short acceleration epoch. The second and final accelerated epoch is of the de Sitter 
type. We do not plot — q(i), Hit), pit), and U(t), see [13] for such plots. We are interested in light propagation and 
unlighted epochs. Fig. [5] illustrates this class of models. The unlighted epoch harbors now the first and second epochs 
of deceleration, as well as the first epoch of accelerated expansion. Clearly, this is a new feature. 




FIG. 4: One transition point universe. This is an example of a universe with one transition point in its evolution. The plots in 
panels (a), (b), (c), and (d) are presented for a typical model with the following parameters: £=0.1, a=50, V( )=l, -E( O )=0.0205, 
A( ) = l, p„=0.333T0~ 4 , and II'(1) =— 5. The plots of all functions presented in the panels (a), (b), (c) and (d), look like the plot 
for the hyperbolic tangent. The plots have no extrema and illustrate the one-fold transition from a universe expanding with 
negative acceleration to a universe with accelerated expansion. The history of this universe is clearly divided into two epochs 
of acceleration/deceleration without distinguished eras inside. This is clearly seen in panel (a) for H. During the first epoch of 
deceleration an unlighted epoch can arise, see panel (b) for which Q = 1 was chosen, and both phase and group velocities are 
less than the speed of light in vacuum. The second epoch can be characterized by parameters which are close to cosmological 
values measured nowadays. In particular, the phase and group velocities during all the second epoch, see panels (c) and (d) 
where Q — —0.35 was chosen, are close to one. 
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FIG. 5: Two transition points universe. This gives an example of a universe evolution with two transition points and one 
extremal point. The plots in panels (a), (b), (c), and (d) are presented for a typical model with the following parameters: 
5=0.1, cr=-0.08, V (0 ) = l, £(o)=0.0205, A (0) = l, /9,=0.333 ■ 10~ 4 , and lT(l)=-15. Panel (a) displays the function H. In panel 
(b) one can see that the plot for the square the refraction index n 2 (t) contains an unlighted epoch of the second type since the 
function n 2 (t), here for Q = —0.28, vanishes and then takes negative values at the end of the first era of the first epoch, in 
contrast to the case, when an unlighted epoch appears at t = 0. In panel (c) the phase velocity is plotted. Other new feature 
is that in the early universe the group velocity, calculated here for Q = 0.2, is close to the vacuum speed of light, and then this 
velocity reaches the same value asymptotically at t — > oo, see panel (d). 



6. Quasiperiodic universe 



The sixth class of models is the class of universes with a large albeit finite number of transition points. It is a 
quasiperiodic universe. This model is intermediate between the model with two transition points and the periodic 
model. 

This class arises when, starting from some transition point, the curve — q(t) remains above the line q=0. This means 



that at a later time the universe's expansion is accelerated. We do not plot — q(t), H(t), pit), and II(t) , see [17j for 
such plots. 

We are interested in light propagation and unlighted epochs. Fig. [7] illustrates this class of models. The behavior 
of the functions H(t), n 2 (t), Vph(i) and V gT (t) is quasiperiodic and the amplitudes of their oscillations decrease 
asymptotically. In addition to the unlighted epoch of the first type, one can see now one point with n = 0, i.e., a 
compressed unlighted epoch of the second type. In between them the square of the refraction index is positive. The 
behavior of the phase velocity is unusual in the early universe since it is greater than the speed of light in vacuum 
during the first era of the first epoch, and then tends to the asymptotic value equal to one in a quasioscillatory manner. 



7. Generic features of the six models: The photon pathlength, the true duration of epochs and the universe lifetime 
Let us define two dimensionless functions, 

i r* 

T ph (t) = - / dt'V ph (t'), (64) 
1 Jo 

and 

r gr (*) = i J dt'V SI (t'). (65) 

The first of them, r p h(i), gives the ratio between two quantities: the length traveled during a time t by an elec- 
tromagnetic wave with phase velocity V p h(t) and the length traveled during the same time t by a photon moving 
with the speed of light in pure vacuum. In the second quantity, r gr (i), the phase velocity is replaced by the group 
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FIG. 6: Three transition points universe. The parameters of the model are the following: Q 



for panels (a) and (b), 



-0.5 for panels (c) and (d), £=0.1, a=l, V( )=l, £ ( o)=0.0205, A (0 )=l, p»=0.333 • 10 , and n'(l)=-5. The universe 
passes through two epochs of deceleration and two epochs of accelerated expansion. The unlighted epoch harbors now the first 
and second epochs of deceleration, as well as the first epoch of accelerated expansion. 
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FIG. 7: Quasiperiodic universe. The universe history is divided into a great number of epochs by a finite number of transition 
points. Panel (a) of H(t) shows clearly this behavior. The plot of n 2 (t), see the panel (b) with Q = 2, demonstrates that in 
addition to an unlighted epoch of the first type, arising at t = 0, an unlighted epoch of the second type compressed into a 
point also appears. This model shows also an unusual behavior of the phase velocity in the early universe, namely, it is greater 
than the speed of light in vacuum during the first era of the first epoch, in contrast to the models described above, and then 
tends to the asymptotic value equal to one in the quasioscillatory regime. The behavior of the group velocity is of a different 
kind, namely, in the early universe it grows monotonically, and then tends asymptotically to one in a quasioscillatory regime, 
see panels (c) and (d) with Q — 0.3. 



velocity. Alternatively, one can consider these two functions as mean values of the phase and group velocities averaged 
over time t, respectively. We have calculated these two functions for the models discussed above. The results are 
presented in the Fig. [3] A general feature of all presented plots is that as t — > oo both these functions tend to one 
limj-^oo r p h(i) = 1 and lim t _ i . 00 T gr (t) = 1. This is not surprising, since when t — > oo the contributions of the periods 
with V p h 7^ 1 and V gr ^ 1, which in themselves are short, become vanishingly small. Thus, for the estimation of the 
total life time of the universe we can assume that, on average, photons propagate with speed of light in pure vacuum. 
On the other hand, when we have to calculate the duration of epochs and eras in the early universe, the behavior of 
T p h and r gr , as well as the estimation of the photon pathlengths, depend on the type and parameters of the model, 
and on whether one uses the phase velocity or the group velocity. 
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FIG. 8: Solid lines illustrate the behavior of F ph (i) and dashed lines illustrate the behavior of r gr (t). r ph (i) gives the ratio 
between the phase length, which the electromagnetic wave runs during a time t, if we take into account its phase velocity, and 
the length traveled by a photon moving with the speed of light in pure vacuum. For F gr (i), the phase velocity is replaced by 
the group velocity. Panel (a) illustrates the periodic model, panel (b) relates to the model with perpetual acceleration, panel 
(c) contains the illustration to the model with one transition point, and panel (d) illustrates the quasiperiodic model. For the 
periodic and quasiperiodic models the plots of the functions F p h(i) and F gr (t) have no intersection points. The other two models 
show the presence of crossing points. Generally, in the late-time universe the quantities T ph (t) and F gr (i) give practically the 
same results, the ratios are close to one. The difference between these quantities is, however, essential in the early universe. 



V. CONCLUSIONS 



We can draw several conclusions. 

1. The standard, concordant, cosmological model deals with at least two epochs, during which electromagnetic 
waves cannot scan the universe's internal structure neither bring information to observers. The first epoch is when 
photons are in local thermodynamic equilibrium with other particles, and the second is when photon scattering by 
charged particles is strong. One can call these two periods of cosmological time as standard unlighted epochs. After 
the last scattering surface, photons become relic photons and turn into a source of information about the universe. 

Now, if one takes into account electromagnetic interactions with the dark sector, i.e., with the dark energy and 
dark matter, one can expect that unlighted epochs of a new type can appear. Here we described one possible 
example, namely, the unlighted epochs produced by the non-minimal coupling of gravitational and electromagnetic 
fields. Since in the framework of the non-minimal three-parameter Einstein-Maxwell model, the curvature coupling 
can be formulated in terms of an effective refraction index n(t), as we did, we can take advantage of the well-known 
classical analogy, namely, in a medium with n 2 < electromagnetic waves do not propagate and their group velocity, 
i.e., energy transfer velocity, has zero value at the boundary of the corresponding zone. Keeping in mind such an 
analogy we studied here, both analytically and numerically, cosmological models admitting unlighted epochs, the 
photon coupling to the spacetime curvature being the key element of the models. 

2. We established a formula for the group velocity, or energy transfer velocity, of an electromagnetic wave non- 
minimally coupled to the gravity held, namely, V gr (t) = ^"/^ , with |V^ r (i)| < 1. Since unlighted epochs with a 
negative effective refraction index squared can start and finish only when n 2 (t^)=0 or n 2 )=oo, the group 
velocity is zero at the unlighted epochs boundary points. At these points the electromagnetic energy transfer stops, 
and thus the condition V gv (t)=0 is the condition to be employed as the criterion for recognizing the beginning and 
the end of the unlighted epochs. 

We have distinguished three types of unlighted epochs: the first one starts at t—0 with negative effective refraction 
index squared and finishes when n 2 —0; the second one starts at t=tm > and finishes at t~t^ 2 ) > with 
n2 (^(i)) =n2 (£(2))=0; the third one is characterized by the fact that at least at one of the two boundaries n 2 has an 
infinite value. The phase velocity V p h{t), another important physical characteristic of an electromagnetic wave, is such 
that V p h(t) = -^j, and so, when n 2 =0 it can take infinite values at the beginning or at the end of unlighted epochs. 
The appearance of unlighted epochs of all three types were illustrated analytically, using a cosmological model with 
a scale factor of the Kohlrausch-type, i.e., a stretched exponential scale factor. 

The unlighted epochs of the first and second types were then described numerically using our Archimedean-type 
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interaction model between dark energy and dark matter. Clearly, the appearance or absence of unlighted epochs 
is connected with the signs and absolute values of the non-minimally coupling parameters qi, q 2 and q^. Thus, 
fingerprints of the unlighted epochs in the cosmic microwave background data could give some constraints on these 
non-minimal coupling parameters. We hope to discuss this problem in another work. 

3. The appearance of unlighted epochs caused by the non- minimal coupling of photons to the gravitational field 
adds new features into the history of the universe as written by electromagnetic fields. Note that if the universe passes 
through an unlighted epoch of the second type, we know for certain that information scanned by electromagnetic waves 
during the preceding epochs is washed out, i.e., such unlighted epochs act as informational laundry. 

In this connection there is a very interesting question: When and why unlighted epochs of the second type can 
appear? Let us note, that quite generally, the behavior of the function n 2 (t) inherits the features of the function H(t); 
to see this compare the panels (a) and (b) in Figs. [2H3 I n this sense, the example of the model with two transition 
points, see Fig. [5J shows explicitly that the appearance of an unlighted epoch of the second type is a consequence 
of large quasioscillations of the function H{t) and thus a consequence of large quasioscillations of the spacetime 
curvature. Such quasioscillations, in their turn, are the result of the dark matter and dark energy energy-momentum 
redistribution, induced by the Archimedean-type coupling fl6l [l7[ . 

4. The time span of different eras in the universe history can be estimated using optical information. From these 
durations one can estimate the distance traveled by the electromagnetic waves through their velocity of propagation. 
But then, one should clarify what is the electromagnetic wave propagation velocity which should be used in the 
interpretation of the observational data. Is it the group velocity, the phase velocity, or the velocity of light in pure 
vacuum? As we have seen these quantities differ in the early universe in an essential way. From Fig. [5] it is seen, 
for instance, that the group velocity and the corresponding traveled distance are less sensitive to the details of the 
dynamics of the universe. 
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Appendix 

As mentioned in subsection III C 11 here we give several examples of specific non-minimal theories by motivating 
possible choices for the non-minimal coupling parameters q\, q 2 and q^. Using the two effective non-minimal coupling 
constants, Q\ and Q 2 given by 

Qi = -2{3 qi +2q 2 +q 3 ) , Q 2 = -2{3 qi +q 2 ) , (66) 

we have shown we can write the square of the refraction index in terms of Hubble function H (t) and acceleration 
parameter —q(t) as 

2 _ l + [Q2~Qiq(t)}H 2 (t) 

H[) l + [Ql-Q2q{t)}HHt)- [b '> 

Of course n 2 (t) depends on Q\ and Q 2 , and thus on qi, q 2 and q 3 . However, independently of this choice, for 
the de Sitter universe, one obtains n 2 (t) = 1. Indeed, since the acceleration parameter is constant and is equal 
to one, —q(t)=l, we obtain that e(t) = — 1+ (Q2+Q1) H 2 (t). Although e(t) 7^ 1 and /j,(t) 7^ 1, one has 

n 2 {t) = e(t)/x(t) = 1 for arbitrary coupling parameters. 

In general, the cosmological models are sensitive to the signs and the values of the parameters Q\ and Q 2 . There is 
an infinite variety for the choice. However, six phenomenological possibilities using geometric analogies and physical 
motivations can be implemented. 
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Vanishing non-minimal susceptibility scalar (1Z — 0) 



When 6gi+3(?2+93 = and thus, from Eq. (|28|). 7Z = 0, one has Q\ + Q 2 — 0. Putting Q = Qi = —Q 2 , one obtains 
in this one-parameter example the following expression, 

2 _ l-QH 2 (t) + m l + QH(t) 
n[t) -l + QH*(t)[l + q (t)]-l-QH { t)- ( ' 



This example has attracted our attention and we have used it in our numerical calculations in Section IIVI Indeed, 
the interest in this choice is that only one function, H , guides the behavior of the effective refraction index, with 
n 2 — > 1, when H —t 0. From a physical point of view this means that asymptotically electromagnetic waves coupled 
non-minimally to the gravitational field propagate with a phase velocity equal to the speed of light in vacuum. 

Gauss-Bonnet-type relation 

When one imposes that the differential equations forming the non-minimal Einstein-Maxwell system are of second 
order (see, e.g., [20|, HI]), then one should use qi+q 2 +q3 = and 2q 1 +q 2 = 0. In this case the non-minimal suscep- 
tibility tensor IZikmn is proportional to the double dual Riemann tensor * R* kmn , i.e., IZikmn = l*R*kmn f° r some 
constant 7 [43j]. Then we have Qi — 0, Q 2 = —2qi, and 

s l-2 qi H 2 (t) , x 

(69) 

If the non-minimal parameter q\ is negative, then during the acceleration epoch (i.e., —q{t) > 0) the square of the 
refraction index is positive, and there are no unlighted epochs. 

First Weyl-type relation 

Assume now that 1Z mn — 0, in which case one can take 3qi + q 2 = and q 2 + 93 = 0. Then, the susceptibility 
tensor is proportional to the Weyl tensor Cikmn, i-e., TZikmn = uCikmn, and Qi — 0, Q 2 = 0. In this case the effective 
permittivity scalars are equal to the unity, e(t) = 1 and /i(t) = 1. Thus n 2 (t) = 1 in any epoch. 

Second Weyl-type relation 

Take 3q\+q 2 = and q% = 0. Then the susceptibility tensor is proportional to the difference between the Riemann 
and Weyl tensors, i.e., IZikmn — Q[Rikmn — Cikmn]- In this example Qi=6qi, Q 2 =0, and thus 

2 l-6q iq (t)H 2 (t) 

71 {t) = TT^hW) ■ (70) 

Now the square of the refraction index is positive during the acceleration epoch (—q(t) > 0) when the non-minimal 
parameter qi is positive. 

Symmetry relation with respect to the left and right dualizations 

If one imposes that the left and right dual tensors coincide *fz lkmn =TZ*' lkmn , then one gets the example in which 
92+93 = 0. In this example Q\=Q 2 , and thus n 2 (t) = 1 for arbitrary epochs, although the dielectric permittivity and 
magnetic permeability 

e{t) = -L = 1 - 2(3 9l + q 2 )H 2 [1 - q(t)} (71) 



differ from unity when 3gi+g2 7^ 0. 
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Drummond-Hathrell-type relation 

An example of a calculation for the three coupling parameters based on one- loop corrections to quantum electrody- 
namics in curved spacetime has been considered by Drummond and Hathrell [21] . In this example the non-minimal 
coupling parameters are connected by the relations q\ = — 5Q, q 2 = 13Q, = —2Q. The parameter Q is positive, 
and constructed by using the fine structure constant a and the Compton wavelength of the electron A e , i.e., Q = 1SQ " n . 
So, Qi = -18Q, Q 2 = 4Q, yielding 

n2 l + 2QH\t)[2 + 9g(t)}_ 
l-2QH 2 (t)[9-2q(t)} 

This example was discussed in [38j . 
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